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Abstract 

We study the lowest lying 7r+ 7r“ resonance R, bound by the long range 
Coulomb potential and destabilized by short range strong interactions 
mediating the dominant decay into two neutral pions. Parametrizing 
the corresponding partial decay width by 

r ( i? ^ 27r° ) = I a 3 y/ ml+ - { A a d m^+ ) ^ 

the quantity A a d depends within QCD & QED on four basic param¬ 
eters : Aqcd , mu+rud , a , rud - mu- 

We are interested here only in the limit where the last two parameters 
vanish 

A a°jj = A a D i ^qcd , mu+md , a = 0 , Wd - = 0 ). 

While isospin invariance implies A a ^ = a — a it is shown 
that beyond the first order expansion in the strong interaction the iden¬ 
tities a = a ^, where a ^ denote the strong interaction scattering 
lengths, do not hold. 


^Contribution to the Int. workshop ’’Hadronic atoms and positronium in the standard 
model” , Dubna, 26-31 May 1998. 

^Work supported in part by the Schweizerischer Nationalfonds. 





1 Introduction, general and specific issues 


The general physics situation addressed here deals with the interplay of a 
long range but weak binding force ( specifically Coulomb attraction within 
7r“ mesonic atoms or pionium ) and strong but short range forces ( repulsive 
and/or attractive ) which however do not themselves give rise to binding 
(specifically destabilising Coulombic bound states, which become resonances 
due to the strong transition 7r“X —> T in particular 7r“7r'’“ —> vr^ vr®). 

In this paper we will exclusively deal with the lowest lying (pionium) 

resonance, denoted R, which in purely nonrelativistic Coulomb spectroscopic 
terms corresponds to the IS state of the composing charged pions. Yet the 
above general situation not only applies also to 7r“p and 7r“d systems, i.e. 
pionic hydrogen and deuterium i , i , but applies as well to (almost) 
classical gravitationally bound systems, e. g. Halley’s comet, subject to the 
short distance solar radiation pressure, which eventually will dissolve the 
bulk of the comet upon successive near approaches. 

For our system R, the (partial) width T { R —> 2Tr^ ) is given by the 
expression 

Qf 

T (R ^ 27r^) = fB - \ T2^r\^ 

8 TT m ^ 

\/ - 4 m 2 ^ (1) 

qf = _ ; fs = \ '■ Bose factor 

2 

T2^r = {2■, ki,k2\T\R]p) ; p = k i + k 2 

In eq. (||) state vectors are labelled by particle content in asymptotic out¬ 
going plane waves, characterized by four momentum vectors k 1^2 and p , 
subject to energy momentum conservation. The relativistic normalization 
of one particle states is adopted throughout. 

\ R ] p ) denotes the pionium resonance state with total energy momentum 
vector p , which is unique modulo redundancies of relative order 
r r/ niR ~ 0.7 10 -9. 

T denotes the momentum conservation reduced T-matrix or transition op¬ 
erator 
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S - 1 


- =T ] { b ; p b \ T \ a ; Pa ) = {2TT)'^6^{pb-Pa) 

i 

{h ] p b \ T \ a ■, Pa) 

( 2 ) 

It follows from the structure of the expression for the decay amplitude in eq. 
i) , that the strong i.e. the short range interactions dominate the transition 
matrix T on one hand, but also modify the structure of the resonance. The 
latter is described by a convenient wave function, best known from nonrela- 
tivistic interactions, represented by potentials in all relevant channels. But 
a relativistic Bethe Salpeter wave function or any variant thereof exhibits 
the same modification in principle. 

The region, where the short range forces modify this wave function (in con¬ 
figuration space) is small compared to the main Coulomb dominated volume 
proportional to the cube of the pionium Bohr radius r^r = 2 j { a m ). 
But this reduced volume is the dominant region from where the decay of the 
resonance takes place. 

To lowest order in the strong interactions the resonance wave function re¬ 
mains unmodified, but if higher orders are included, this modification be¬ 
comes important. Beyond this lowest order no obvious relation exists, even 
though higher order isospin asymmetries and electromagnetic corrections are 
neglected, which expresses the matrix element in eq. in terms of exclu¬ 
sively the strong scattering amplitude of constituent pions on one hand and 
the purely Coulombic wave function for the resonance on the other hand. 
Precisely such a relation has been derived by Deser, Goldberger, Baumann 
and Thirring Q for pionic atoms. Since its validity would encompass the 
much more general interplay of long range weak binding and short range 
strong but non binding forces, it can easily be falsified in potential mod¬ 
els. The fact, that such models may not be applicable to pionium or more 
generally to pionic atoms is in this respect irrelevant. 

In section 2 we will derive such a relation involving beyond the abovemen- 
tioned quantities the dependence of the strong interaction scattering lengths 
on an appropriate coupling parameter ( A ) for strong interactions 

d 

A a 5) = A - A a ° 

d A 

A a ° = ( a - a -^=2 ) ( Aqcd, mu + m^, a = 0, = 0 ) 

(3) 
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In eq. (^) the quantities A a ^ , a ^ refer to the scattering length in the 
limit specified. 

Estimates of the lifetime of pionium are presented in section 3. 


2 Resonance decay amplitude in the Lippmann- 
Schwinger framework 

The amplitude T 2 ^ _r in eq. (|^ is proportional to the reduced amplitude 
in relative space coordinates as adapted to the center of mass frame 


R = K + \ H^ \ R; rel") ; k = - k 2 ) 


( 4 ) 




2 , rel ; k + 


In eq. (Q) 1 denotes the interaction Hamiltonian in the Schrodinger repre¬ 

sentation. The state | R ; rel ) is normalized according to the scalar prod¬ 
uct in relative coordinate space. Because of the open decay channel the 
resonance is an eigenstate in the continuous spectrum of the total relative 
Hamiltonian, including all electromagnetic interactions 


H rel \ R ] rel )= m R \ R ; rel ) ; H rel = H 0 + H i H 

( 5 ) 

Furthermore ip in eq. (^) denotes the outgoing scattering wave function 
obeying the Lippmann-Schwinger equation 



1 

E - ie - H 


( 6 ) 


E = m R = 2 -s/ k + rn ^ ; { x \ ip j: ) = exp { i k x ) 

In eq. @ V' g denotes a plane wave. 

A fully relativistic description of relative coordinate space involves the use 
of relative time and the corresponding Bethe Salpeter equation . 

Finally the kinematic constant K in eq. (Q) is determined from the equivalent 
expression for the resonance width in eq. (|^) 

r ( fi ^ 21“ ) = fa j ( 2 ir) S {E - m 


k 

R ) - 

( 27r ) 3 
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From eq. (|^ we determine the constant K in eq. 

Qf m R 


r ( i? ^ 27r« ) = /b 


4 TT 

K = \/2 ( m B ) 3/2 


(^ip i \ H I \ R ; rel^ 


( 8 ) 


Associated 2 ' —> 2 scattering amplitudes 

The resonance decay amplitude in eqs. © and @ is associated with the 
scattering amplitude for the 2^ —> 2 pion reaction 


vr^ TT^ 


TT ' vr 


da {2' 2 ) 


k 1 k 


= fs 


Pi P 2 


d n 


2^2 


9/ 


di 


( 9 ) 


64 TT ^ s 


2^2 


/ = (2 7r°;A;i,/c2|r|7r+7r ; p i , p 2 ) 


All kinemaric quantities in eq. (^) refer to the center of mass system. 

The invariant amplitude T 2 ^ 2 ' depends on the standard Lorentz invari¬ 
ants 

s = {p 1 + p 2)'^ ; t = {p 1 - k i) ; u = {p I - k2)‘^ 


ki + k2=Pi+P2 = {E,d) 


9/ 


2 


Qi 



2 


( 10 ) 


Analogous to the relative coordinate space decay amplitude 

rel ^ in eq. (|^ the corresponding 2 ' ^ 2 scattering 

amplitude is 


T 2 ^ 2 ' — 


Hi 1 

^p) 

; k = 

i ( 1 - k 2 ) 


■ft) = 

2 7r° , rel ; 





(11) 

II 

- P 2 ) ; 

1 V' p 

> = 

vr"'' 7r“ , 

rel ; p ) 
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I V' P ) in eq. (0) refers to the plane wave asymptotic vr^ vr state. 
The expression for the cross section in eq. (^) becomes 


d 3 k 


da [2' ^ 2) = fB {2 t:) 5 {E - E') 


(2vr) 




^ lab 


E = 2 J k + m I ^ ; E' = 2 Jp'^ + m 


( 12 ) 


2 p E 


^ lab — 


E — 2 m 


Thus we obtain the expression for the differential cross section equivalent to 
eq. (I) 


da {2' 2) 


d O 


fs 


df 


Qi 


16 vr 


E ^ — 2 m 




Parametrizing the diferential cross section in the form 


da {2' 2) 


d fl 


= Ib 


<lf 


Qi 


F 


and comparing eqs. (|^ and (|14D it follows 
1 


F = 


vr E 


T 


2^2 


m 


4 vr 




2 p 


1 + 




m 


(13) 


(14) 


( 15 ) 
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From eq. (p!5|) we deduce the kinematical constant K 


K' 


2 m TT ^ E 


1 + 




2 p 


2 


m 


2 

TT + 


E = yfs 


(16) 


The nonrelativistic limit involves the shift of the energy E —> E — 2 m ^ 
to zero value at vr"*" 7r“ threshold. 


3 Relations between decay and scattering ampli¬ 
tudes 

We extrapolate the rate formula from above threshold to the resonance 
position with an arbitrary incident intensity I 

I = V lab Q ; Q = \ { E ) 

d a {2' ^ 2) 

T {E Q ^ 2t^^) = dn - / 

dn 


= Ib 


Qf E 


4 vr 


\ Eil'ipp) X 


(17) 


to E —> m R , i.e. below vr'*' 7r“ threshold. It follows, that the resonance is 
described equivalently by the appropriate choice of the amplitude X ^ X R 
from eq. (|^ 

ip i \ H I \ R ] rel = x R E ^ m r p t \ H i \ -ip 


— 4 m ^ 


(18) 


P 


< 0 


We define the extrapolated scattering length a r below 7r+ tt threshold in 


accordance with eq. QlSD 


4 TT 


lim E ^ m R ip i \ H I \ ijj = - a R V2 


m 


(19) 


aR = aR(”27r°”^7r+7r ) 
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The combinatorial factor \/2 included in the definition of a in eq. (p^ ) 
accounts for the two identical pions in an algebraic way. The quotes in 
” 2 vr ^ ” refer to the bose symmetrized and normalized state. 


Substituting a in eq. 
the resonance 


(17) we obtain upon extrapolating to the mass of 


T ( i? ^ 27r° ) = 


4:TT qf m R 


a RX R 


( 20 ) 


m 


-IT + 


It is rewarding in view of sequential approximations to use as unit of density 
the Coulomb density | X c I ^ for the vr'*' 7r“ system 


m 


1/2 


X c = \ a 


; C R = 


8 TT 


X R 


X c 


( 21 ) 


The expression for the resonance width in eq. (2C) becomes 

m R 


T { R ^ 2Tr^ ) = Qf 


m TT + o. R i R 


2 m 


( 22 ) 


a D = a R C R 

In eq. (^2|) the quantity a d denotes the decay equivalent scattering length, 
to be distinguished from a r which is indeed the scattering length, extrap¬ 
olated to the resonance mass and for the reaction ” 2 <— vr'*' 7r“. 

Up to this point all equations were exact, to be evaluated in QCD , QED 
theory, restricted to two light quark flavors with masses m u , iR d- 
Hence a ^ is a function of 4 basic parameters 


a D = a D { ^ QCD ,mu + md, a, md-ruu) (23) 
We are only interested here in the limit, where the last two of the four 


parameters in eq. (^) a and m d — Rt- u tend to zero, whereby only the 
lowest order contributions to the resonance width T { R —> 27r*^ ) are 

retained : 


a , m d - m u 0 ; ^ m ^ 

r(7?^27r°)^r°(7?^27r°) ; m r ^ 2 m ^ ^ 

an ^a% = aD{^ qcd, mu + Rid, a = 0,md-rn 


0 ) 

(24) 
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The momentum variables p , k also tend to zero in the above limit according 
to eq. (|I^. 

Before going over to the systematic approximation deffined in eq. 
represent the resonance width T ( i? —> 2-71^ ) in the form 

T {R ^ 2Tr^) = r 0 ( ^ ^ 27r° ) ( 1 + J ) 


we 


6 = 6 { a , m d - m u) = 0(a, (m d - m u ) ^ ) 


(25) 


In eq. (^) only the relevant expansion parameters of the correction factor 
6 are exhibited explicitly. 

The dependence of 5 dominated by the first order a correction has been 
discussed in this workshop by A. Rusetsky and H. Sazdjan i > 0 and 
amounts numerically to 


0.06 


(26) 


27r0 


enhancing in lowest nontrivial order the limiting width T *^ ( R 
to which exclusive attention is directed in the following. 

The limiting situation : a , m d— 'mu 0 

From eqs. ( p^ and ( p^ we deduce the following expression for the limiting 
width T ° ( R —> 27r° ) 

n I 2 


{R 


0 _ 


D 


= a 


2TT^) = a 


R R 


- m\ 


m 


D 


(27) 


In the limiting situation u-d isospin is an exact symmetry, which implies the 
general decomposition into 1 = 0,2 channels remembering the channel 


definition in eq. 

(HD 





n 0 - 

^ D , R — 

n 0 

^ D , R 

to 

o 

T 

^TT ) 



TT ■’“vr “ 

) = 

O 

II 

nr 
V 3 

\I = 2) 

(28) 

” 2 vr 0 

”) = 

vT u = 0) 

+ 

^ V 3 

\I = 2) 


It follows from eq. (|^) relinquishing the subscript R 
length 

on o , the scattering 

O 

T 

: 0° 

v/2 

— A a ° ; 

3 

n 0 - 

a D — 

%/2 

A a D 

3 

(29) 

A o ° = o 

— a 

^=2 ; Aa% = 

n -^=0 - 
a 

a ^=2 
a 
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In the following we drop also the superscript on the quantities 
A a ° —> A a and A a 5) —> A a d • 

Substituting eq. in eq. ( |^ ) we obtain 

r°(i?^27r°) = |a^ ^ - m\ | m ^ A a d | ^ (30) 


For the purpose of getting orders of magnitude into perspective we define 
the reference width, substituting the quantity A a for A a in eq. (|30|) 


r , ( i? ^ 27r0 ) = 


= ^ a 3 

9 “ 


- rn\ | m ^ A a 


ref 


= ^ a ^ J rn\^ - m\ 


(31) 


r „ ( ^ 27 r 0 ) = r 


ref 


m . A a 


0.25 


Substituting the reference width eq. (pH) becomes 


m 


A 


r ( i? ^ 27 r 0 ) = r ref 


0.25 


A a D 


A a 


T ref = ( 5.397 ) 10 ^ - m\ = 0.1917 eV 

T ref = I / T ref = ( 3.434 ) 10 “I® SeC 


(32) 


We focus on the limiting behaviour of the quantities (^(pi\Hi\R] rel ^ 

in eq. (^) in comparison with the scattering amplitude (p t \ H i \ 
in eq. (H) and their limiting relation in eq. 


if i \ H I \ R ] rel^ 




4 vr 

I - X C ^ a D 


m 


TT + 


(33) 


4 TT 

I - X c A a 

m n ^ 


Follwoing S. Deser et al. |Q] we consider the purely Coulombic bound state 
of TT + vr “ , which we shall denote \ ip q ) as unperturbed state together 
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with the limiting perturbing Hamiltonian 


Hi H 

and the eigenstate(s) | 'h ) of the full Hamiltonian 

H = He + ; H c = H 0 + Vc 

From the two equations 


H \^) = E \^) 

; He 1 

V' c ) = 

Ec 

1 C ) 

A E = E — E c : 

II 

o 

2 m ^ 

1 

4 

a ^ m 7r 


(34) 

(35) 


(36) 


the limiting relation for the energy shift A E oi the resonance follows 


{'ll; C \ H I ^ 


A E = 


( "0 c I ^ ) 

In the limit we are considering the energy shift becomes 

A E = {^c \ H f - 

\ Ec - H 


(37) 


^ str. 


V’ c 


V 2 
A C 


Re T ( TT ■*“ TT 


vr ' TT 


4 m 


m 


= —a 


m ^ [2 a ^ ~ + a'^“2 


(38) 

It is instructive to follow S. Deser et al. |Q] and extend the relations in eq. 
to the (complex) elastic tt + vr “ scattering amplitude 


1 


P 


Ec - H 


1 


E c — H — i e 


A E - i'y / 2 = - xc 
X = T a{R ^ 2 tt^) 


T {tT TT \ TT TT ) 


4 m 


( 39 ) 
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The quantity 7 = T a { R —> 27r° ) in eqs. ( |^) and (|3^) is the limiting 
width of the state | '0 C ) as it decays, to all orders in the strong interaction, 
to 2 vr ° . 


This corresponds to the substitution in eq. 


R ; rel ) 

I'lp c) 


l^ip t \ H 1 \ R ] rel'^ 

^ (ipi\H,\TPc) 

(40) 

A a D 

—> A a 



However the state | R ; rel ) as it evolves from without external pertur¬ 
bation to a time t = 0 say, is not | V' c ) j but rather the full incoming 
Schrodinger state \ f ~ ;C ) , adiabatically evolving from \ tp c ) 


1 


ip ;C ) = \ Ip c ) + 


E PI — H — i e 


jj str. 


V' c 


(41) 


with the energy E ji including the energy shift A i? in eq. (^) . 
In the limit we are considering it follows 


1 ; 

C) 


Xc\v~]tt^tt-) 


( + 1 

Hi 1 

R ; rel ) 

XC 





((/?+; 2 TT \ H 1 \ If - 

] TT ^ TT ) 

; = p 

= 0 

H - 

- Hi 



(42) 

From the structure of limiting amplitudes in eq. (42) we infer , dropping 
the superscript in the following 


1 


'Ip ] ”2 TT 


55 Q O55 


H 


I - GoHi 

Go = {Ethr — Ho + i^) 


Ip ; TT TT 


(43) 


We recall that ip in eq. (43) refers to plane wave states with vanishing 
momentum at threshold. 


The argument for the limiting substitutions in eq. (|42D is like this ; in the 
matrix element {(p~^\Hi\R; rel ) the action of H i —> H restricts 
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the configurations composing both states \ ^ ) and | R ; rel ) to normal 

strong interaction relative distances d ~ 1 ~ 4 fm. 

Those configurations are in the sense of the limit considered insensitive to 
the two key parameters governing the resonance at its determining distance, 
i.e. its Bohr radius 
a B = 2/(m7 r+«) ~ 400 fm 
As a consequence also the key mass square difference 
— ~ ( 35.51 MeV ) ^ 

plays no significant role. Hence the dominant configurations relevant for the 
transition amplitude within the state | R ; rel ) are the same as those in 
which R is (almost) bound , i.e. for a / 0 but m = m . Oi course 
the resonance R is always decaying into two (and more) photons. 

In the above situation the (almost) bound state R is by no means described 
by the Coulomb wave function, in particular at the distances within d rel ■ 

It follows combining eqs. (^3|) and (|^ 


— A a 


D 


3 m 


8 TT 


ijj ; 2 TT ^ \ H I 


I - Go Hi 


'0 ; vr TT 


— A a 


3 m 


8 TT 


V’ ; 2 TT ° I H 1 


1 - Go Hi 


ijj ; TT ^ TT 


(44) 


shows that decay and scattering amplitudes are not the same. 


Eq. 

A a B and A a according to eq. (|44D can not be related to each other 
without detailed knowledge of H i . 

Let us introduce the coupling strength A - always remaining in the 2 fla¬ 
vor SU2 symmetric QCD limit, with fixed m ^ through the 

substitution 


Ho , Go 

- Ho,Go 


Hi 

A Hi 

(45) 

A a D , A a 

-> Aa£)(A),Aa(A) 
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Then it follows from eq. (^) 


Aa£)(A) = A 


d A 


A a ( A 


(46) 


In the limit considered the quantities A a d , A a depend within QCD on 
the two basic parameters A 2 qcd and mu = m d ■ 

An equivalent set is / tt , the pion decay constant , and m ^ 


= m 


It follows from the relation in eq. (^) that the pion mass is to be held 
constant, whereas at least in lowest two orders of chiral perturbation theory 
IP] the variation of A is equivalent to a variation of / 


d 


d 


d A 


d Q 


(47) 


m 


Q = k 


; k arbitrary fixed constant 


f: 


Estimates of pionium lifetime 

to estimate A a d and the lifetime of pionium, 
m n 
are 


We use eqs. (^) and 
with the shorthand notation m 


TT + 


The lowest order (tree level) values 


7 m 


m TT a 


/ = 0 _ 
( 1 ) 


m TT a [ 1 ) 2 = 


32vr/2 

2 ml 

32 vr/; 


= 0.1562 


= -0.0446 


(48) 


To this end we list from Ecker et al. 0 the contributions through two loop 
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order to both scattering lengths a ^ ^ and a ^ ^ 


order 

m a ^ 

T — 2 

m a 

m A 

1 

0.16 

-0.045 

0.205 

2 

0.04 

0.003 

0.037 

3(1) 

0.017 

0.0007 

0.0163 

3 {II) 

0.006 

-0.0023 

0.0083 

total (/) 

0.217 

-0.0413 

0.258 

total (II) 

0.206 

-0.0443 

0.250 


We base our estimate on the one loop contribution to a ^ ® which is of the 

form 


/ = 0 _ 


( 2 ) 


= k I Q ^ L ; L = — log Q 


(50) 


The numerical value of the logarithm in the one loop contribution to a ^ ^ ^ 
is quite accurately 4 as in eq. (^) , when all nonlogarithmic terms of order 
Q ^ are absorbed into the argument of the logarithm. The constant k in eq. 
(^) can then be chosen such that q = exp —4 according to eq. (^) . 

It then follows 


d 


d A 


/ = o _ M _ r-M ^/ = o ^ 3^1 = 0 


( 2 ) 


= ( 1 - L ) 


( 2 ) 


I a ^2)= = 0.03 (51) 


Neglecting all other corrections from higher orders and from the 1 = 2 chan¬ 
nel we obtain as our estimate 


A a 


D 


A a ~ 0.28 ; A a ~ 0.25 
With the ratio A a d / ^ cl given in eq. 
r 0 ( i? ^ 27r0 ) = 1.25 T ref 

= 0.240 eV 


- ~ 1.12 

A a 

we obtain from eq. 

A ^ 2 

m jr + A a 


(52) 


0.25 


m 7r + A a 


(53) 


0.25 
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or equivalently for the lifetime, using the abbreviation F ( i? —> 27r^ ) 

r 0 


0.25 


To = 1 / r 0 = ( 2.74 ) 10 sec 


m TT . A a 


(54) 


Finally we apply the radiative corrections as estimated by A.Rusetsky and 
H.Sazdjan @ , Q according to eqs. (^) and (^) and find for the resonance 
width the estimate 


r ( R ^ 27r0 ) 


1.33 r 


ref 


= 0.255 eV 

or equivalently for the lifetime 
r ( R ^ 27r° ) = ( 2.58 ) 10 sec 


m TT A a 
0.25 

m TT A a 
0.25 


1 + <5 
1.06 
1 + 
1.06 


(55) 


0.25 


1.06 


m TT + A a 


1 + 5 


(56) 

If we omit the correction proportional to{AaD/Aa)‘^ the corresponding 
estimate for the lifetime would be 


0.25 


1.06 


T {R ^ 2tt^) = { 3.24 ) 10 -15 


sec 


m 


A 


1 + 5 


(57) 


The results on the above lifetime to date are from L. Nemenov et al. 01 
T { R ^ 2^“ ) = ( 2.9 ^ ) 10 sec (58) 


In conclusion we are looking forward in suspense to the measurement or 
better the analytic deduction of the lifiteme of pionium from the study of 
the breakup reaction in targets with appropriately chosen thickness by the 
DIRAC collaboration |T^ . 
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